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Cretan matrices are orthogonal matrices with elements < 1. These may 
have application in forming some new materials. There is a search for Cretan 
matrices, especially with high determinant, for all orders. These have been 
found by both mathematical and computational methods. 

This paper highlights the differences between theoretical and computational 
solutions to finding Cretan matrices. 

It has been shown that the incidence matrix of a symmetric balanced 
incomplete block design can be used to form Cretan(u; 2) matrices. We give 
families of Cretan matrices constructed using Hadamard related difference 
sets. 

Keywords: Hadamard matrices; orthogonal matrices; Cretan matrices; symmetric 
balanced incomplete block designs (SBIBD); difference sets; 05B20. 

1 Introduction 

Cretan matrices were first discussed, per se, during a conference in Crete in 2014 by 
N. A. Balonin, M. B. Sergeev and colleagues of the Saint Petersburg State University 
of Aerospace Instrumentation, 67, B. Morskaia St., 190000, St. Petersburg, Russian 
Federation but were well known using Heritage names, [1, 2, 8, 9]. This paper follows 
closely the joint work of N. A. Balonin, Jennifer Seberry and M. B. Sergeev [3, 4, 7]. 

We highlight the difference between mathematical solutions to the problem and com¬ 
putational solutions to the problem where round off errors can be crucial. We stress that 
in real world applications in engineering and building it it impossible to cut physical 
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materials precisely. Moreover changing climate conditions make some error tolerance 
absolutley necessary. 

An application in image processing led us to search for T-variable orthogonal matrices, 
S, with maximal or high determinant. We use x\, X 2 , •••, x T as variables. When the 
variables are replaced by entries/values/numbers having modulus < 1 and at least one 1 
per row and column they are called Cretan matrices. The entries/values/numbers can be 
negative in the mathematical theory but only non-negative in practical applications. 

Hence the aim of our study is to find T-variable orthogonal matrices which yield Cretan 
matrices which have maximum or high determinant and as few variables as possible by 
using both mathematical and computer aided strategies. 

We extensively use the article by Seberry [14]) for definitions and the mathematical 
approach. 

Symmetric balanced incomplete block designs or (v, k, A)-configurations or 
SBIBD(v,k, A) are of considerable use and interest to image processing (compression, 
masking) and to statisticians undertaking medical or agricultural research. We use the 
usual SBIBD convention that v > 2k and k > 2A. 

We see from the La Jolla Repository of difference sets [13] that there exist (v, k, A) 
difference sets for v = At + 1, 4 t, 4t - 1, 4t - 2 which can be used to make circulant 
SBIBD(v, k, A). 

In this and future papers we use some names, definitions, notations differently to how 
we they have been have in the past [2]. This we hope, will cause less confusion, bring our 
nomenclature closer to common usage, conform for mathematical purists and clarify the 
similarities and differences between some matrices. We have chosen to use the word level, 
instead of value for the entries of a Cretan matrix, to conform to earlier writings [2, 8, 9]. 

1.1 Preliminary Definitions 

Although it is not the definition used by purists we use orthogonal matrix as follows. 
Definition 1. A orthogonal matrix, S, has real entries and satisfies 

S J S = SS J = L0l n 

where I n is the n x n identity matrix, and u, the weight , is a constant real number. For 
computationally discovered cases we will say it is orthogonal to (say) five decimal places. 

Definition 2. S = ( Sij ) of order n will be called a T-variable orthogonal matrix , with vari¬ 
ables x\,X 2 ,"', x T when it is orthogonal, satisfying SS 1 = wl and for which Y?j=i{ s ij ) 2 = u, 
oj a real constant, for all i and X!/=i s ]i s ki = 0 for each distinct pair of distinct rows j and 
k. A similar condition holds for the columns of S. We write this as S = CM(ji',t ; •••). 

Notation 1. In this paper we only study 2-variable orthogonal matrices of order n, S 
or S(v\ 2; •••), written with the variables x and y. 

When the variables are replaced by real numbers with modulus < 1, (these may 
be negative numbers), the resultant matrix is orthogonal and called a Cretan matrix 
CM(n; 2; w; ■■•). The original 2-variable matrix and the resultant orthogonal matrix S 
are used to denote one-the-other. 
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Definition 3 (Cretan). A Cretan matrix, S = (s tJ ), of order n, with r levels, written 
as Cretan(n;r) or CM(n;r',uj), is a r-variable orthogonal matrix, which has had the 
variables replaced by real numbers with modulus < 1, and CM(n;r; ca) and has a least 
one 1 in each row and column. A CM(n;r), S, satisfies the orthogonality equation 


S T S = SS J =ul n , (1) 

here uj, called the weight, is a real constant. 

A Cretan matrix is either precisely orthogonal or orthogonal to (say) 5 decimal 
places. □ 

One repeating question is to find how close a computational approach comes to a 
precise approach and to find out the real world meaning of two or more solutions which 
arise in the mathematical approach. Certainly the extra solutions can not have maximal 
determinant but what do they tell us? 

Definition 4 (Levels). The entries of the r-variable orthogonal matrix of order n, S, 
are called variables. When the variables are replaced by real numbers/values/entries with 
modulus < 1, and there is at least one 1 in each row and column we have a Cretan(n\T\u ) 
or CM(u;t-,uj ) which has r levels. 

That is the number of different real numbers, counting plus entries separately from 
minus entries, is the number of levels of S. 

The level x = 1 is pre-defined for all Cretan matrices. □ 

We emphasize: a Cretan(n;r) or CM(n;r) or S, have r levels, they are made from 
t- variable orthogonal matrices by replacing the variables by appropriate real numbers 
with moduli < 1, where at least one entry in each row and column is 1. 

Notation 2. S, a r-variable matrix with variables x\, X 2 , •••, x T may be used to form 
a Cretan(n;r; uj-,ji,j 2 ,---,j T ',^i,i 2 ,-"i^T'ideterminant): where n is the order, r is the 
number of distinct variables or levels (counting -x separately from x); .)i,i 2 ,'",.)r is the 
number of occurrences of the variables xi,X 2 ,—,x t if they occur the same number of 
times in each row and column, however as this mostly does not happen these values 
are just omitted; 4 ,^ 2,"’]4 is the total number of each variable in the whole matrix 
Cretan(n, r) or CM(n,r), and the determinant. 

After variables have been replaced by feasible entries/values/numbers S, Cretan(n,r) 
or CM(n,r ) or S, are used, loosely to denote one-the-other. 

Cretan matrices may be used to find some real matrices with entries < 1 which have 
with maximal or high determinant. In this paper r = 2. 

For 2-level Cretan matrices we will denote the levels/values by x,y where 0 < \y\ < x = 1. 
We also use the notations Cretan(v), Cretan(v)-SBIBD and Cretan-SBIBD for Cretan 
matrices of 2-levels and order v constructed using SBIBDs. 

A Cretan(n]T]UJ ) matrix, or CM(n;r;uj) has r levels or values for its entries. 

For r-level Cretan matrices we will denote the levels/values as {l\,l 2 , ■ ■ ■ ,C-}- More 
generally we can have notation such as 
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1. CM (order ), 

2. CM (order] r = number of levels), 

3. CM (order, t = number of levels; u =weight), 

A. CM (order; r = number of levels; u =weight; determinant), 

5. CM(order; t = number of levels; u ^weight ; levels = ji,j 2 , ■ ■ ■ occurrences of 
levels in whole matrix = £i, £ 2 ,,l T ), 

6. CM (order; t = number of levels; ui = weight ; levels = ji, J 2 , • ■ •, j T )', occurrences of 
levels in whole matrix = £i,£ 2 , ■ ■ ■ ,£ T ; determinant), 

7. 

according to the parameters of current importance. The definition of Cretan is not that 
each variable occurs some number of times per row and column but each variable occurs 
£ 1 , £ 2 , ■ ■ ■, £t times in the whole matrix. So we can have CM(re; r; co; levels; £i,£ 2 , ■ ■ ■, £ T \ 
determinant). Q 

1.2 Notation Transitions 

In transiting from one mother tongue to another (Russian to English and English to 
Russian) and from previous to newer usage, some words reoccur: we need a shorthand. 
To simplify references we note: 


Heritage Usage 

Cretan Matrix 

References 

Fermat 

CM(4t + l) 

[2, 4] 

Hadamard 

CM( At) 

[3, 9, 4] 

Mersenne 

CM (At - 1) 

[1, 2, 9, 4, 5, 6] 

Euler 

CM(At - 2) 

[1, 8, 9], 


Usage 

Dual Usage 

CM(4t) 

CM(4t-3) 

Core of CM (At + 1) 
Core of CM (At - 2) 


The odd order Hadamard matrices and the Mersenne matrices referred to in [2] and 
[1] are SCM(4t+l) and SCM(4t-l). 

1.3 The Theoretical Precise Case 

We now define our important concepts the orthogonality equation, the radius equation(s), 
the characteristic equation(s) and the weight of our matrices. 
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Definition 5 (Orthogonality equation, radius equation(s), characteristic equa¬ 
tion^), weight). Consider the matrix S = ( Sij ) comprising the variables x and y. 

The matrix orthogonality equation 


S T S = SS J = ul n 


yields two types of equations: the n equations which arise from taking the inner product 
of each row/column with itself (which leads to the diagonal elements of coI n being ui) are 
called radius equation(s), g{x , y ) = oj, and the n 2 - n equations, f(x, y) = 0 , which arise 
from taking inner products of distinct rows of S (which leads to the zero off diagonal 
elements of t oI n ) are called characteristic equation(s). The orthogonality equation is 
Ylj=i s 1j = w. u is called the weight of S. □ 

Example 1. Exact example 

We consider the 2-variable S matrix given by 


x 


S = 


V 

y 

y 

.y 


y 

x 

y 

y 

y 


y y 
y y 

x y 

y x 
y y 


y 

y 

y 

y 

X. 


By definition, in order to become an orthogonal matrix, it must satisfy the radius and 
characteristic equations 

x’ 2 +4y 2 = cj, 2xy + 3y 2 = 0. 

To make a Cretan(5;2;y) we force x = 1, (since we require that at least one entry 
per row/column is 1), and the characteristic equation gives y = — Hence cu = 3|. The 
determinant is (^)a = 20.286. Thus we have an S = Cretan(5 ; 2 ; y; 20,5; 20.286). □ 


1.4 The Algorithmic Computations Case 

Hybrid-precise-computational. 

Example 2. Consider the 9-variable S matrix given by 


w 1 

£l 

£2 

£3 

£4 

£4 

UJ 2 

£l 

£2 

£3 

£3 

£4 

UJ 3 

£l 

£2 

£2 

£3 

£4 

CJ 4 

£l 

. £1 

£2 

£3 

£2 

w 5 . 


where cu* « 1 - Si, i = 1, ■•■5 and e j & e « 0, j =1, 2, 3, 4, to (say) 5 decimal places and Si 
and £j are very small numbers. Then the orthogonality equation S a S ^ = cu/ 5 , the radius 
equations and the characteristic equations can be specified exactly and solved to (say) 5 
decimal places. S a is a Cretan(5;9;cu) hybrid-precise-computational matrix. 
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Now consider the case for CM(5 ; 2 ; u>) 


c 


S b = 


b 
-a 
-a 


b 


b -a -a b' 
c b -a -a 
b c b -a 
-a b c b 
-a -a b c 


with radius equation c 2 + 2 b 2 + 2a 2 = ui and characteristic equation 2 be + a 2 - 2 ab = 0. 
These can be solved for n = 5 and 


a = 1.000000, , b = 381965, c +0.309017 


giving 


SrSl 


'2.387286 

0.000001 

0.000002 

0.000002 

. 0.000001 


0.000001 

2.387286 

0.000001 

0.000002 

0.000002 


0.000002 

0.000001 

2.387286 

0.000001 

0.000002 


0.000002 0 . 000001 ' 
0.000002 0.000002 
0.000001 0.000002 
2.387286 0.000001 
0.000001 2.387286. 


Since all the equations can be specified exactly they can be solved exactly (or with very 
small errors) or to as many decimal places required. It is a hybrid-precise-computational 
Cretan(5; 3; 2.387286) matrix. 


Computational example. 


Example 3. The following computationally discovered matrix S c has maximal deter¬ 
minant, 3.3611175556, and is a Cretan(5; 3; 3||) matrix, order 5 when a=l.000000, 
b=0.500002, and c=0.333340. In fact it has a precise solution a=l, b=| and |. 


a 

-a 

-b 

-a 

-c 

b 

a 

-a 

c 

-a 

a 

a 

c 

-b 

a 

a 

-c 

a 

a 

-b 

c 

-b 

-a 

a 

a. 


The main questions that arise with finding Cretan matrices from computer calculations 
are 

• Can we find the equations for the exact solution? 

• Can we find solutions with very small errors? 


2 Preliminary Definitions and Results: SBIBD 

Definition 6 (Incidence Matrix). For the purposes of this paper we will consider 
an SBIBD(v, k, A), B , to be a v x v matrix, with entries 0 and 1, k ones per row and 
column, and the inner product of distinct pairs of rows and/or columns to be A. This is 
called the incidence matrix of the SBIBD. For these matrices X(v - 1) = k(k - 1). 
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We note that for every SBIBD(v, k, A) there is a complementary SBIBD(v, v-k,v- 
2k + A). One can be made from the other by interchanging the 0’s of one with the l’s 
of the other. The usual use SBIBD convention that v > 2k and k > 2A is followed. For 
examples see [14]. 

In this work we will only use orthogonal to refer to matrices comprising real elements 
with modulus < 1, where at least one entry in each row and column must be one. 
Hadamard matrices and weighing matrices are the best known of these matrices. We 
refer to [4, 15, 12, 17, 14] for definitions. 

2.1 Mathematical Foundations for the 2-Variable Orthogonal Construction 

Let S' be a 2-variable matrix of order n: S will be written with variables x\,X2,--,x s 
S = S(n;r;cu; ji,j s ;^i,^ 2 )"')^s;determinant) where n is the order, r is the number 
of distinct variables (counting -x separately from x)\ 

j\. j ‘2 , ■■■■j.s is the number of occurrences of the variables xi,X 2 ,---,x s if they occur the 
same number of times in each row and column, (however as this mostly does not happen 
these values are just omitted); 

^ 1 ,^ 2 ) is the total number of each variable in the whole matrix S, and the 
determinant. The original 2-variable matrix and the resultant orthogonal matrix S after 
the variables have been replaced by feasible entries/values/numbers are used to denote 
one-the-other. 

In all these Hadamard related cases (v = At - 1) (but not necessarily in all cases) the 
2-variable orthogonal matrix with higher determinant comes from the SBIBD(At- 1,2 1, t ) 
while the SBIBD^At - 1, 2t — 1, t - 1) gives a 2-variable orthogonal matrix with smaller 
determinant. These examples are given because they may give circulant SBIBD when 
other matrices do not necessarily do so. 


3 Mathematical vs Computer Aided Constructions 


In all the mathematical constructions the search is for a precise solution. However 
in engineering and building it is not possible to manufacture to such precision. Both 
the measurements of the physical materials used and the computer solutions will have 
tolerance built into their outcomes. Computer aided solutions necessarily come with 
round-off errors in the calculations. The aim then is to find solutions which are of 
practical use. 

To a mathematician an orthogonal matrix X will satisfy XX 1 = ul. However for real 
world use the matrix 


x 


S = 


y 

y 

y 

.y 


y 

X 

y 

y 

y 


y y y 
y y y 
x y y 
y x y 

y y x. 
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written as 



'w 1 

£1 

£2 

£3 

£ 4 ' 



£4 

cu 2 

£l 

£2 

£3 


= 

£3 

£4 

cu 3 

£l 

£2 

5 


£2 

£3 

£4 

W 4 

£1 



. £1 

£2 

£3 

£2 

w 5 _ 


where cu* « 1 -<£*,* = 1, ••■5 and e.j * 

e » 

■o* 

0 

= 1 , 

2, 3, 

4, to (say) 5 decimal places and 5, 


and ej are very small numbers is a practical solution. The the orthogonality equation 
S a S[ = iu/ 5 , the radius equations and the characteristic equations can be specified exactly 
and solved to (say) 5 decimal places but is computationally very time-consuming. S a is a 
Cretan(5;9;o;j) hybrid-precise-computational matrix. 

3.1 Two Computational Maximal Cretan Matrices: The Balonin-Mironovski 
Matrices, (BMC) 

We illustrate in Figure 1 two famous examples [2] the Balonin-Mironovski Ay and the 
Balonin-Mironovski Ag both with maximal determinant. We use the notation BMC(order; 
number of levels = r; weight = u ; levels; ; determinant) ( or sometimes BMC(order; 
number of levels = r; weight = uj ; levels; total occurrences of each of the r levels; 
determinant). These matrices are of special interest as they are close to a singular point 
in a computer program investigating metal combinations. 

The Balonin-Mironovski-Cretan(7;5;5.0777;30,6,3,4,6) uses a first row and column 
{d b b b a a a} around a three block circulant core, see http://mathscinet.ru/catalogue/ 
definitions/ for details. 

This is quite different from the Balonin-Mironovski-Cretan(9;4;6.4308;40,16,24,1) which 
uses a first row and column {-d bbbbbbbb} around a circulant core circ(a -a c c a 
c -a -a). A 7 and Ag have maximal determinant for Cretan(7) and Cretan(9) matrices 
respectively. 

We do not have a theoretical model for these two cases as shown in Figure 1. 


4 Precise Orthogonal Matrices from SBIBD 

We now use SBIBD s to construct 2-variable orthogonal matrices from SBIBD(x,y) s. 
We always, in making 2-variable orthogonal SBIBD(x,y ) from an SBIBD, change the 
ones of the SBIBD into x and the zeros of the SBIBD into y. 

We use the Main SBIBD 2-variable construction theorem from Seberry [14] for 
our mathematical construction. Note there are two parts one from the orginal design and 
the other from it’s complement. The solutions may appear the same but arise differently. 

Theorem 1. [Two-level Cretan Matrices from SBIBD:I] Let S be made from an 
SBIBD(v, k , A), B, by replacing the 1 ’s with x and the 0 ’s with y. Then S is a Cretan(v; 
2; w) or CM(v; 2; u) where 

u>i = kx 2 + (v - k)y 2 (2) 
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(a) BM A7: BMC(7;5;5.0777;30,6,3,4,6) (b) BM A9: BMC(9;4;6.4308;40,16,24,1) 

Figure 1: Balonin-Mironovski matrices A7 and A9 with Maximal Determinant 


and the characteristic equation is 

Xx 2 + 2 (k - X)xy + (v - 2 k + A )y 2 = 0 . (3) 


The determinant is ut) . 

We see this leads to precise numerical answers. Using the complementary SBIBD(v , v- 
k, v - 2k + A) we have 

Theorem 2. [Two-level Cretan Matrices from SBIBD.TI] Let S be made from 
an SBIBD(v, v - k, v - 2k + A), B, by replacing the 1 ’s with x and the 0 ’s with y. Then 
S is a Cretan(v; 2; 002 ) or CM(v; 2; 0 J 2 ) where 

0 J 2 = (v- k)x 2 + ky 2 (4) 

and the characteristic equation is 

(v-2k + X)x 2 + 2 (k - X)xy + Xy 2 = 0 . (5) 


y_ 

The determinant is uj f . 


Corollary 1. [The Cretan 2-level Matrices from SBIBD Theorem], or [ The 
Cretan-SB IB D(v; 2) Theorem] Whenever there exists an SBIBD(v,k, A) there exist 
two Cretan(order;T;w,y,x ; determinant), or 2-level Cretan-SBIBD(v; 2), or S or CM, 
as follows, 


1. 


One from the Cretan(v\2;kx 2 + (v-k)y 2 ; ^ k ^ 2 fc+C ~~^ 
the SBIBD(v,k,X ) and 


1; determinant); made from 
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2. one of the Cretan{v\ 2 ; kx 2 + (v - k)y 2 ; ^ ^ ’ 1 5 determinant), made from 

the SBIBD(v,k ,\) or 

the Creta/nfv: 2; ( v-k)x 2 +ky 2 ; p determinant), made from SBIBD(v, 

v - k, v - 2k + A). :JEJ 

4.1 Mathematics for Some Hadamard Matrix Related Constructions 

There are three obvious Hadamard related constructions (but these are by no means all): 
those using SBIBD(4t- 1, 2t-l, t- 1), those using the Menon difference sets and those 
using the twin prime difference sets. We illustrate using the first. 

Corollary 2 (From Hadamard Matrices). Suppose there exists an Hadamard matrix 
of order 4 1, then there exists an SBIBD(4t — 1, 2t-l, t - 1). 

Hence for an S(v = 4 1- 1;2), satisfying Equations ( 2 ) and (3) variables we have 
a Cretan(4t-l;2;uj\; y, x; determinant) from the solution for the 2-variable orthogonal 
matrix, by setting x = 1 in 

x, y = -1^1 = “—~~ - w i = 2t+(2t-l )?/ 2 and det(S') = {2tx 2 +{2t-\)y 2 )^r . 

For S = S(v; 2), satisfying Equations (4) and (5) we have a Cr et an (4t-1; 2; W 2 / y > x; 
determinant) matrix from the solution for the 2-variable orthogonal matrix, by setting 
x = 1. In particular for the SBIBD(4t - l,2t,t) we have 

1 , y= ~j~ Y~’ 1^1 = 1 ~ 1; ^2 = (2t-l) + 2ty 2 and det(S') = (2t-l + 2ty 2 )^ 1 . 

Explicitly for SBIBD(v, k, A) = SBIBD{ 4f-1,2t — 1, t-1): v = 4t — 1, k = 2t — 1, A = t-1, 
v - k = 2t, k- X = t, v - 2k + A = t, and we have the Cretan matrices: 

• Cretan (4 1 - 1; 2; 2t - 1 + 2 ty 2 - =^, 1 ; (2t- 1 + 2 ty 2 )^) ; 

• Cretan ^4t - 1; 2; 2t - 1 + 2 t?/ 2 ; , 1 ; ( 2 1 - 1 + 2 ty 2 )^~ j . □ 

In all these Cretan - Hadamard cases (but not in all cases) the Balotin-Sergeev- 
Cretan(v) matrix with higher determinant comes from the SBIBD(4t - 1 , 2 t,t) while 
the SBIBD(4t - 1, 2t — 1, t - 1) gives a Cretan(At - 1) matrix with smaller determinant. 
These examples are given because they may give circulant SBIBD when other matrices 
do not necessarily do so. 
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(a) the principal solution (b) the complementary solution 


Figure 2: Orthogonal matrices for order 7: Balonin-Sergeev Family 

To construct solutions of the second type (where the weighing matrix is an attractor), 
we choose the SBIBD, B, to be that with the smaller number of l’s per row and column, 
i.e., v > 2k. 

Corollary 3 (Menon Sets and Regular Hadamard Matrices). Suppose there 
exists a regular Hadamard matrix of order 4m 2 , then there exists an SBIBD(4m 2 , 2m? - 
m, m 2 -m). Hence we have two-level Cretan matrix, S, satisfying Equations (1) and (3) 
for y = ^?jx. The principal solution, from the SBIBD(4m 2 , 2m? - m,m 2 - m) is well 
known as a regular Hadamard matrix with 

x = 1, y = 1, u\ = 4 m? has det(S) = (4m 2 )"” 1 . 

The second solution from the SBIBD(4m 2 ,2m? + m, m? + m) gives 

*= 1 . y-M*’ det(S) = ((^) 2ro ” 

for a two-level Cretan matrix with smaller determinant. □ 

Use is made of these matrices a small study in [7]. 
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(a) H (n = 16) (b) S (n = 16) 

Figure 3: Regular Hadamard and 2-variable orthogonal matrices for order 16 


Example 4. Consider p = 3 which gives an SBIBD( 13,4,1). To generate this matrix we 
use circ(- 6 , a, - 6 , a, a, - 6 ,- 6 , -b, a, - b , - 6 , -b,-b), with characteristic equation a 2 - Gab + 
6b 2 = 0. 

There are two solutions, see Fig. 2a where b = 3+ ^ , and Fig. 26 where 6 = 3 ~^ //3 , so 
we have 

a = l, 6 =—j_ = 3 *n/3 cj = 4o 2 + 96 2 = 7± — , = 9.5981 or 4.4019 

3tV3 6 2 

13 

for the required two-level orthogonal matrices. They have det(5) = ^7 ± 2 = 

2.4221 x 10 6 or 1.5264 x 10 4 . □ 

We note the appearance of some of the matrices constructed in more than one family 
that has been studied: this family with two-level Balonin-Sergeev-Mersenne [9] and the 
Cretan-SBIBD-Singer matrices discussed above. □ 


5 Conclusions 

In all the mathematical constructions the search is for a precise solution. However in 
Engineering and Building it is not possible to manufacture to such precision. Both 
the measurements of the physical materials used and the computer solutions will have 
tolerance built into their outcomes. Computer aided solutions necessarily come with 
round-off errors in the calculations. The aim then is to find solutions which are of 
practical use. 

The many questions arise with finding Cretan matrices from computer calculations. 
Besides those questions mentioned above, where exact equations can be found to match 
computational solutions, can we use structured matrices, such as those, for example, 
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which are circulant, two circulant, have a border, have two borders, to yield valuable 
information towards finding precise solutions, without any errors or with very small 
errors? 

We see from the La Jolla Repository of difference sets [13] that there exist (v, k, A) 
cyclic difference sets and hence SBIBD(v , k, A) for v = At + 1, At. At - 1, At - 2, where t 
is integer, which can be used to make circulant SBIBD(v,k, A). We recall that the two 
CM(At - 1) where v = 3 (mod 4) 2-level matrices arise, one from the SBIBD and the 
other from its complement. This result is not necessarily so for v in other congruence 
classes. 

We note the existence of some CM(At - 1) 2-level matrices with the same parameters 
from the Cretan(At - 1) - Singer family and the Balonin-Sergeev (Mersenne family) [9], 
which are defined for orders At - 1, where t is integer [8] . We have not considered the 
equivalence or other structural properties of CM matrices with the same parameters. All 
other useful references to this question may be found in [4], 

Matrices of the Cretan(At + 1) - Singer family also have orders belonging to the set of 
numbers At + 1, t odd: these are different from the three-level matrices of Balonin-Sergeev 
(Fermat) family [8, 9] with orders At + 1, t is 1 or even. The latter exist for orders v, 
where v - 1 is order of a regular Hadamard matrix, described above via Menon sets. 

Orders At + 1, t is odd, are Cretan(At + 1) - SBIBD matrices; their order may be 

2 

neither a Fermat number (2 + 1 = 3, 2 2 + 1 = 4 + 1, 2 2- + 1 = 16 + 1, 2 22 + 1 = 256 + 1, ...) 
nor a Fermat type number ( 2 even + 1). The first CM( 37; 3) is discussed in [7]. It uses 
regular Hadamard matrices as a core and have the same (as ordinary Hadamard matrices) 
level functions. We call them Cretan(At + 1) - SBIBD matrices and will consider them 
further in our future work. 

The twin prime power difference sets allow us to have circulant Cretan(At-l)-SBIBD 
in orders pq which were not peviously known, pq = 35, 143, 323, 

The main conclusion (about alternating matrices) follows: orders Ak + 1, k odd, 
belong to alternating two- and three-level matrices of the Cretan-SBIBD-Singer and 
Balonin-Sergeev-Fermat families. The sets of orders known for Cretan-SBIBD-Singer 
and Balonin-Sergeev-Mersenne families have orders in common, they (and hyper-plane 
SBIBD based matrices) are similar to each-other and distinct from the three level 
Balonin-Sergeev-Fermat family. 

The unexpected main conclusion is that the two Cretan(v, k ) we find by this method 
can either 

1. both arise from the SBIBD(v,k, A); 

2. both arise from the SBIBD(v, v - k, v - 2k + A); 

3. one arises from the SBIBD(v,k, A) and the other arises from the SBIBD(v,v- 
k, v - 2k + A); 

Cretan matrices are a very new area of study. They have a cornucopia of research 
lines open: what is the minimum number of variables that can be used; what are the 
determinants that can be found for Cretan(n; r) matrices; why do the congruence classes 
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of the orders make such a difference to the proliferation of Cretan matrices for a given 
order; find the Cretan matrix with maximum and minimum determinant for a given 
order; can one be found with fewer levels? how can computational constructions help or 
be helped find optimal of near optimal solutions to problems? 

We conjecture that will give unusual conditions. □ 
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